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ABSTRACT 



As part of a far-infrared (FIR) spectral scan with HerschelfPACS, we present the first detection of the hydroxyl radical (OH) towards 
' the Orion Bar photodissociation region (PDR). Five OH (X ^H; v=0) rotational A-doublets involving energy levels out to £„/A;~51 1 K 

a] have been detected (at -65, -79, -84, -119 and -163yum). The total intensity of the OH lines is 2^(OH)=i5xlO"'' ergs"' cm"- sr"'. 

The observed emission of rotationally excited OH lines is extended and correlates well with the high-y CO and CH+ J=3-2 line 
(~| ■ emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models 

including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (T/^^ldO- 
I I 220 K) and dense (nn^lO^^^ cmT^) gas with source-averaged OH column densities of >10'^cm"^. High density and temperature 

O ■ photochemical models predict such enhanced OH columns at low depths {Av<V) and small spatial scales (-10" cm), where OH 

formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH 
C/5 ' emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps 

I or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much 

higher than those expected in equally warm regions without enhanced FUV radiation fields. 

Key words, astrochemistry — infrared: ISM — ISM: abundances — ISM: molecules 
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\^ 1 . Introduction Besides, OH is the product of H2O photodissociation, the main 

pO ' . . , destruction route of water vapour in the gas unshielded against 

^^.[In quiescent regions iiradiated by cosmic- or X-rays, the oxy- radiation. Observations of OH in specific environments 

in gen chemistry is initiated by the charge transfer from H+ and H+ ^j^^^ constrain different chemical routes of the oxygen chemistry. 

to atomic oxygen, forming O and OH . In warm environments 
^""^ it can also start with the reaction of atomic oxygen with H2(v=0) 



Previous observations with KAO and ISO have demonstrated 



^ . to form OH. This endothermic reaction (by -0.08 eV) possesses *at OH is a powerful tracer of the warm neutral gas in shocked 

> ■ an activation barrier of a few thousand K and high gas tempera- 8^^^ fr""^ protostellar outflows and super n ova reninants to ex - 

irl tures (>400 K) are needed to produce significant OH abundances tragalactic nuclei (e.^., 'Storev et alJUMIt | Melnick et al. || 1987t 

X - (e.g., in shocked gas). In molecular clouds exposed to strong Gonzalez-Alfonsoetal. 2004). Unforfiinately, the poor angular 

b ■ far-ultraviolet (FUV) radiation fields, the so-called PDRs, the resolution (>1 ) and sensitivity of the above telescopes prevented 

- ^ ■■ gas is heated to relatively high temperatures and there ai-e also ^om resolving the OH emission from intersteUar PDRs. 

high abundances of FUV-pu mped vibrationally excited m olec- Owing to its proximity and nearly edge-on orientation, the 

ulai- hydrogen H5(v=l,2...) (iHollenbach & TielensI [19971) . The interface region between the Orion Molecular Cloud 1 (OMCl) 

internal energy available in H* can be used to overcome the and the Hii region illuminated by the Trapezium cluster, the 

OeP) + H2(v=0) reaction baiTier rsee lAgundez et1Dl20T()l and O^ion Bar, is the prototypical wai-m PDR (with a FUV radia- 

references therein). Although not well constrained observation- tion field at the ionization front of ^:^2.5xl0 times the mean 

ally, enhanced OH abundances are expected in warm PDRs. mtersteUar field in Draine units; Marconi et al. 1998). The most 

OH is a key intemiediary molecule in the FUV-illuminated commonly accepted scenario is that an extended gas component, 

gas because further reaction of OH with H,, C+, O, N or S+ mean gas densities of 10' cm ; , causes the chemical 

leads to the formation of H2O, CO+, O2, NO or SO+ respectively, stratification seen in the PDR d Hogerheijde et al.| 11995D. Most 

of the low-y molecular line emission arises in this extended 

* Herschel is an ESA space observatory with science instruments "interclump" medium ( Tielens & Hollenbach 1985; Simon et aD 

provided by European-led Principal Investigator consortia and with im- 119971 Ivan der Wiel et al . 2009; Habart et al.ll2010l) . In addition, 

portant participation from NASA. another component of higher density clumps was introduced to 
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Fig.l. OH (X ^n; v=0) rotational lines detected with Herschel/PACS towards the q'2ooo^5''35"'21.9', 5200o^-5''25'06.7" position 
where the higher excitation OH line peak is observed (see Figure|2]i. The red dotted lines show the expected wavelength for each A- 
doublet (see Figure lBTl in the appendix|B]for a complete OH rotational energy diagram). Small intensity asymmetries are observed 
in most OH A-doublets. Transition upper level energies and A/j spontaneous transition probabilities are indicated. 



fit the observed H2, hig h-7 CO, CO ^ and other high density 
and temperat ure tracers (Burton et al. 1990; Parmar et al. 199lt 
Storzeret all 119951: Ivan der Werf et al...l996; Young Owl et all 
Owing their small filling factor this clumpy structure 
would allow FUV radiation to permeate the region. The pres- 
ence of dense clumps is still controversial. 

In this letter we present initial results from a spectral scan of 
the Orion Bar taken with the PACS instrument (iPoglitsch et alJ 
I2OIO ) on board Herschel (Pilbratt et al. 20ldi) as part of the 
"HEXOS" key programme (Ber gin et al.l i2010). PACS observa- 
tions of OH lines towards young stellar objects have recently 
been reported by Wampfler et al. (I l2010b . Here we present the 
first detection of OH towards this prototypical PDR. 

2. Observations and data reduction 

PACS observations were carried out on 7 and 8 September 
2010 and consist of two spectral scans in Nyquist sampling 
wavelength range spectroscopy mode. The PACS spectrometer 
uses photoconductor detectors and provides 25 spectra over a 
47"x47" field-of-view resolved in 5x5 spatial pixels ("spax- 
els"), each with a size of ~9.4"x9.4" on the sky. The measured 
width of the spectrometer point-spread function (PSF) is rela- 
tively constant at /l<125jum but it increases above the spaxel 
size for longer wavelengths. The resolving power varies between 
^/A^~1000 (Rl grating order) and -5000 (B3A). The central 
spaxel was centred at Q'2000: 5''35'"20.61', (J2000: -5"25'14.0" tai-- 
get position. Observations were carried out in the "chop-nodded" 
mode with the largest chopper throw of 6 arcmin. Nominal ca- 
pacitances (0.14pF) were used. The integration time was 3.2 h 
for the 1342204117 observation (B2B and Rl) and 2.7 h for 
the 13422041 18 observation (B3A and Rl). The data were pro- 
cessed with HIPE using a pipeline upgraded with a spectral 
flatfield algorithm that reduces the spectral fringing seen in the 
Nyquist-sampled wavelength spectra of bright sources. Figure[T] 
shows the resulting OH lines towards the OH emission peak and 
Figure [3] shows the intensities measured in each of the 25 spax- 
els for several lines of OH, CO, CH+, H2O and [N 11]. In order to 
better sample the PSF and obtain accurate line intensities to be 
reproduced with our models, we fit the OH emission averaged 
over several adjacent spaxels in Section|4](see also appendixlAb. 

3. Results 

Of all the observed OH lines only the ground-state lines at 
~119jum show widespread bright emission at all positions. 
The ~119//m lines mainly arise from the background OMCl 
complex (the same applies to most ground-state lines of other 




'^JSOOO 

Fig. 2. PACS rotationally excited OH ^113/2 7=7/2^5/2 lines at 
~84jum (£'„/A:=291 K) overlaid on the distribution of the CO 
7=6-5 peak brightness temperature (colour i mage) observed 
with the CSO telescope at -11" resolution ( iLis et al.iri998h . 
White contours show the brightest regions of H2 v=l-0 5(1) 
emission (iWalmslev et al.li2000l) . Lo wer intensity H^ extende d 
emission is present in the entire field (Ivan der Werf et all 19961) . 
Violet stars shows the position of the H' ^CN 7=1-0 clumps 
deeper inside the Bar (iLis & Schilkel 12003). Note the decrease 
of OH line intensity with distance from the ionization front. 

species). Figure |2] shows the spatial distribution of the rotation- 
ally excited OH ^03/2 7=7/2^5/2 lines at -84 yum (£„//t=29 1 K) 
superimposed over the CO 7=6-5 peak brightness tempera- 
ture (colour image from Lis et al 1 119981) and over the bright- 
est Ht v=1-0 S{\) lin e emission regions (white contours from 
IWalmslev et al.l I2OOOI) . The emission from the other OH A- 
doublets at -79 and -163 ;um (see Figure IBT2] in appendix iBll 
displays a similar spatial distribution that follows the "bar" mor- 
phology peaking near the Hj v=l-0 5(1) bright emission region 
and then decreases with distance from the ionization front (note 
that Hj shows lower level extended emission with small-scal e 
structure in the entire observed field; Ivan der Werf etal1ll996h . 
The excited OH spatial distribution, however, does not follow 
the CO 7=6-5 emission maxima, which approximately trace the 
gas temperature in the extended "interclump" component. 
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Figure [T] shows the detected OH A-doublets (at ~65, ~79, 
~84, ~1 19 and ~163 yum) towards the position where the higher 
excitation OH lines peak. The total intensity of the observed FIR 
lines is 2/(OH)^5xlO ergs ' cm ^ sr All OH doublets ap- 
pear in emission, with intensity asymmetry ratios up to 40% (one 
line of the A-doublet is brighter than the other). 

Note that the upper energy level of the ^113/2 7=9/2— » 7/2 
transition at ~65 //m lies at E„/A:~5 11 K. The critical densities 
(«cr) of the observed OH transitions are high, ricr^lO^ cm"^. For 
much lower gas densities, and in the presence of strong FIR radi- 
ation fields, most lines would have be en observed in absorption, 
espec ially those in the ^113/2 ladder (Goicoec hea & Cemicharol 
Hence, the observed OH lines must arise in an widespread 
component of warm and dense gas. 

Although our PACS observations do not provide a fully sam- 
pled map, the line emission observed in the 25 spaxels can be 
used to carry out a first-order analysis on the spatial correlation 
of different lines (neglecting perfect PSF sampling, line opac- 
ity and excitation effects). Except for the OH ground-state lines 
at ~119/vm (that come from the background OMCl cloud), we 
find that the rotationally excited OH lines correlate well with the 
high-y CO and CH^ emission but, as expected, they do not corre- 
late with the ionized gas emission. Figure[3](^ovver panel) com- 
pares the observed OH ~84.597//m line intensities with those 
of CO 7=21-20 (£„//t~1276K), CH+ 7=3-2 (EJk~240K) and 
[Nil] 121.891 yum (all observed with a similar PSF) and also with 
the CO 7=15-14 (£„/fc~663K), H2O 3o3-2i2 (£„//t~163K) and 
OH ~163.397yL(m lines (upper panel). This simple analysis sug- 
gests that the excited OH, high-7 CO and CH^ 7=3-2 lines arise 
from the same gas component. It also shows that the emission 
from different excited OH lines is well correlated, while the OH 
and H2O emission is not (within the PSF sampling caveats). 

4. OH column density determination 

Determining the OH level populations is no trivial excitation 
problem. In addition to relatively strong asymmetries in the 
collisional rate coefficient^ between each A-doubling compo- 
nent (Offer & Van Dishoeck 1992), radiative and opacity ef- 
fects (pumping by the ambient IR radiation field and line trap- 
ping) can play a significant role if the gas density is much 
lower than «cr- Here we use a nonlocal and non-LTE code that 
treats both the OH line and continuum radiative transfer (see 
appendix in Goicoechea et al. 2006). The continuum measured 
by PACS and SPIRE in the region (H. Arab et al. in prep.) 
can be approximately reproduced by a modified blackbody with 
a colour temperature of ~55 K and a opacity dependence of 
-0.05(100//!)' (see ' Lis et al.]|1998h . Our calculations include 
thermal, turbulent, and opacity line broadening with a turbu- 
lent velocity dispersion of cr=1.0k ms~' and FWHM=2.355 xcr 
(see e.g., the linewidths measured bv^ Hogerheiide et al.ll995l) . A 
grid of single-component models for different A^(OH), gas tem- 
peratures, densities and beam filling factors were run. The best 
fit model was found by minimizing the ";^f^-value" (see the ap- 
pendixlAlfor its definition). 

' We used collision rate coefficie nts of OH with para- and ortho - 
H2 from Ofi'er & van Dishoeck ^'1992') and Offer et al. i ll994h . 
Strong differences in the intensity of each OH A-doublet component 
due to asymmetries in the colhsional rates between OH and para-Hi 
were predicted {e.g., /(84.597)//(84.420)>/(119.441)//(119.234)>l). 
Asymmetries are significantly reduced when collisions with ortho-Wi 
are included {i.e., in the warm gas) and when FIR radiative pumping 
plays a role. We assume that the H2 ortho-to-para ratio is thermalized 
to the gas temperature, e.g., ~ 1.6 at 100 K and ~2.9 at 200 K. 
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CO 15-14 (173,631) ^H^O 3„3-2i2 (174,626) ^OH 163,397 
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^CO 21-20 (124,193) -'CH+ 3-2 (119,852) ^[NII] 12 1 ,898 

Fig. 3. Line intensity spatial correlations between the OH ^113/2 
7=7/2"— » 5/2^ line at 84.597 fim and lines from other species. 
Intensities are in units of 10"^ ergcm"^ s"' sr"'. 



From the excitation models we conclude that the high 
/(65.279)//(84.596):^0.3 and /(65.132)//(84.420)^0.5 intensity 
ratios in the ^113/2 ladder can only be reproduced if the gas is 
dense, at least n//=n(H)-i-2n(H2) of a few 10^ cm"^. In addition, 
no model is able to produce even a crude fit to the data if one as- 
sumes the average molecular gas temperature (Tt^85 K) in the 
lower density "interclump" medium (Hogerheijde et al. 1995]). 
In the ^ni/2 ladder, the intensity of the ~163 yum lines is sensitive 
to the FIR radiation field in the region through the absorption of 
FIR dust continuum photons in the OH ~34 and ~53 yum cross- 
ladder transitions (^03/2-^01/2 7=3/2^5/2 and 7=3/2^3/2 re- 
spectively). However, the ~163yL(m lines in the Orion Bar are 
not particularly strong and the OH ~34 and ~53 pm lines are 
not present in the ISO spectra, thus FIR pumping does not dom- 
inate the OH excitation. All in all, the best model is found for 
a source of high density (nH<10^ cm"-') and warm gas tem- 
peratures (T|t=160-220K). This temperature lies in between the 
~600K derived in the H2 emitting regi ons (lAl lers et al. 2005) 
and the -150 K derived from NH3 lines (lBatrla"& Wilson.,20031) 
near the ridge of de nse and cooler H'^CN clumps (Ta:^50K; 
ILis feSchilkd 12003 ). In our simple model, the best fit is ob- 
tained for a source of small filling factor (7/ ^10%) with a source- 
averaged OH column density of >10'^ cm"^ (see appendix lAb. 

5. OH chemistry and line emission origin 

We used the Meudon PDR code ihe Petit et alJ 120061: 
iGoicoechea & Le Bourlotl2007h to estimate the OH column den- 
sity in a slab of gas at different densities (nn from 5x10^* to 
10^ cm"^). The adopted FUV radiation field, ;^'=10'*, roughly 
corresponds to the attenuation of the FUV field at the ionization 
front by a column density of A^//-10^' cm"^ in a nn-lO* cm^^ 
medium. This attenuation is equivalent to a spatial length of 
~10''cm, consistent with the observed decrease of excited OH 
emission with projected distance from the ionization front (see 
Figure |2]i. Given the high gas temperature, FUV field and mod- 
erate grain temperature (Ty,.~70-100K) in the regions traced by 
FIR OH, CO and CH^ lines, we neglected molecule freeze- 
out and ice desorption (which are important deeper inside at 
Av>3; Hollenbach et al. 2009). In our models OH is a sur- 
face tracer that reaches its peak abundance at Av<l (Figure |4]i 
where OH formation is driven by the endothermic reaction 
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0(V) + H2 ^ OH + H, slightly en hanced by the OCP) + 
reaction (included in our models; see lAgundez et aTll2010l) . Gas 
temperatures around ~1 000-500 K are predicted near the slab 
surface at Av=0.01 and around ~100K at Av>l. In these H/H2 
transition layers where the OH abundance peaks, the electron 
density is still high (<[CVH]n/f) and hydrogen is not fully 
molecular, with /(H2)=2n(H2)/[n(H)+2n(H2)]^0.5. In general, 
the higher the gas temperature where enough H2 has formed, the 
higher the predicted OH abundance. 

In the Av<l layers, OH destruction is dominated by pho- 
todissociation (OH + /iv — » O + H) and to a lesser extent, by 
reactions of OH with H2 to form H2O (only when the gas tem- 
perature and density are very high). Water vapour photodissoci- 
ation (H2O + hv —> OH + H) in the surface layers limits the H2O 
formation and leads to OH/H2O abundance ratios (>1), much 
higher than those expected in equally warm regions without en- 
hanced FUV radiation fields {e.g. in C-shocks). The lack of ap- 
parent correlation between the excited OH and H2O 3o3-2 n lines 
(see Figure |3]l and the absence of high excitation H2O lines in 
the PACS spectra (only weak H2O 22i-2i2, 2i2-loi and 3o3-2i2 
lines are clearly detected) suggests that the bulk of OH and H2O 
column densities arise from different depths. 

As the temperature decreases inwards, the gas-phase produc- 
tion of OH also decreases. The spatial correlation between ex- 
cited OH, CH^ 7=3-2 and high-7 CO lines is a good signature 
of their common origin in the warm gas at low Ay depths. 

Our PDR models predict OH column densities in the range 
-lO'^cm"^ to ~10'^cm-^ at Av<l for gas densities between 
n//=5xlO'* and 10^ cm"^ respectively. Hence, even if we take 
into account possible inclination effects, high density and tem- 
perature models produce OH columns closer to the values de- 
rived in Section 4 just from excitation considerations (note that a 
precise determination of the gas density would require knowing 
the collisional rate coefficients of OH with H atoms and elec- 
trons). The OH abundance in these dense surface layers is of the 
order of ^lO *' with respect to H nuclei. However, optical depths 
of Av<l correspond to spatial scales of only ~10^^ cm (i.e., 
much smaller than the H'^^CN clumps detected bv' Lis & SchiUcd 
I2OO3L deeper inside the cloud), but we detect extended OH emis- 
sion over ~10'^ cm scales. Therefore, we have to conclude that 
the observed OH emission arises from a small filling factor en- 
semble of unresolved structures that are exposed to FUV radi- 
ation (overdense clumps or filaments). Note that owing to the 
lower grain temperature compared to the gas, the expected FIR 
continuum emission from these clumps will still be below the 
continuum levels observed by Herschel/PACS. 

The minimum size of the dense clumps is ~10'^cm (from 
OH photochemical models) with a maximun size of ~10'^cm 
(from the inferred beam dilution factor). Both correspond to 
<0.2"-2" at the distance of Orion. As an example, H2 photoe- 
vaporating clumps of ~10' ^ cm size h ave been unambiguously 
resolved towards S106 PDR lNoel eFa l. (2005). However, higher 
angular resolution observations (e.g., with 8-lOm telescopes) are 
needed to resolve smaller H2 clumps from the H2 interclump 
emission in the Orion Bar. 

If the observed FIR OH line emission does not arise from 
such a high density gas component, a different non-thermal 
excitation mechanism able to populate the OH ^113/2 7=9/2 
and 7/2 levels would be needed. Two alternative scenarios 
can be explored, at least qualitatively. First, OH molecules 
produced by H2O photodissociation are expected to form 
mostly in the ground electronic and vibrational state but in 
unusually high energy 7>70/2 levels (a few thousands K!; 
Ivan Harrevelt & van Heme3l2000l) . Nevertheless, they will cas- 
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Fig. 4. Gas-phase PDR models of a FUV-illuminated slab 
(Y=10^) of gas with different densities: uh-IQ' , 10^, 10^ and 
5x10"* cm-^. OH column densities are shown as continuous 
curves (left axis) while OH abundances, M(OH)/n/f, are shown 
as dashed curves (right axes). 

cade down radiatively to lower energy rotational ladders 
extremely rapidly. Although very excited suprathermal OH 
7^70/2-15/2 lines have been reported towards the HH211 out- 
flow a nd were interpreted as H2O photodissociation (iTappe et al.l 
'2008"), we did not find any of them in the SWS or in the IRS 
spectra of the Orion Bar (taken and processed by us from the 
ISO and Spitzer basic calibrated data archives). Besides, even if 
photodissociation is the main H2O destruction mechanism in the 
Ai/<1 warm layers, this is not the main production pathway of 
OH (the 0(^P) + H2 reaction dominates). 

Second, experiments and quantum calculations suggest 
that reaction of 0(^7") atoms with H2(v=l) can pro- 
duce significant OH in the v=l vibrationally excited state, 
with OH(v=l)/OH(v=0) population rat ios >1 for moder - 
ate kinetic energies (>0.05-0.1eV; e.g.. iBalakrishnanI l2004 . 
Complementarily, absorption of near- and mid-IR photons 
(from the continuum or from bright overlapping Hj and ionic 
lines) can also pump OH to the v=l state. In both cases, subse- 
quent de-excitation through the v= 1 -0 rotation-vibration band at 
~2.80/vm would populate the OH(v=0) rotationally excited lev- 
els that we observe with PACS. However, the OH v=l-0 band at 
~2.80 jum is not present in the ISO/S WS spectra. Even assuming 
that the 0(^P)+H2(v=l) reaction only forms OH(v=l), signifi- 
cant gas column densities at very hot temperatures (T;(.~2000K) 
will be needed to match the observations if hh-IQi^^^ cmT^ . 
Studying the OH vibrationally-pumping mechanism quantita- 
tively is beyond the scope of this work but these chemical and 
pumping effects could contribute to the excitation of FIR OH 
lines in lower density gas. 

Different scenarios for the origi n and nature of photoevapo- 
rating clumps have been proposed (i Gorti & Hollenbach. .2002b . 
but without a more precise determination of their sizes and den- 
sities it is hard to conclude on any of them. Subarcsec resolution 
observations of OH gas-phase products, e.g., direct observation 
of CO^ or SO^ clumps with ALMA, will help us to assay the 
clumpy nature of the Orion Bar in the near future. 

Acknowledgements. PACS has been developed by a consortium of insti- 
tutes led by MPE (Germany) and including UVIE (Austria); KU Leuven, 
CSL, IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAF- 
IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); lAC (Spain). We thank Darek Lis 
and Malcolm Walmsley for providing us their CO 7=6-5 and 1-0 S(l) maps. 
We also thank Emilie Habart, Bart Vandenbussch and Pien'e Royer for useful 
discussions and the referee for his/her constractive report. We thank the Spanish 



4 



Goicoechea et al.: OH emission towards the Orion BarPDR 



MICE^ for funding support through grants AYA2006- 14876, AYA2009-07304 
and CSD2009-00038. JRG is supported by a Ram on y Cajal rese arch contract. 
We acknowledge the use of the LAMDA data base ISchoier et aiT20051 . 



References 

Agundez, M., Goicoechea, J. R., Cemicharo, J., et al. 2010, ApJ, 713, 662. 
Alters, K. N., Jaft'e, D. T., Lacy, J. H. et al. 2005, ApJ, 630, 368. 
Balakrishnan. N. 2004, JChPh. 121, 6346. 
Batrla, W. & Wilson, T. L. 2003. A&A, 408, 231. 
Bergin, E. A., Phillips, T. G., Comito, C. et al. 2010, A&A, 521, L20. 
Burton. M.G., HoUenbach, D.J. & Tielens, A.G.G.M. 1990, ApJ, 365, 620. 
Goicoechea, J.R. & Cemicharo, J. 2002, ApJ, 576, L77. 
Goicoechea, J.R. et al. 2006, A&A, 456, 565. 
Goicoechea, J. R. & Le Bourlot, J. 2007, A&A, 467, 1. 
Gonzalez-Alfonso, E. et al. 2004, ApJ, 613, 247. 
Gorti, U. & Hollenbach, D. 2002, ApJ, 573, 215 
Habart, E., Dartois, E., Abergel, A. et al. 2010, A&A, 518, LI 16. 
Hogerheijde, M.R., Jansen, D.J., van Dishoeck, E.E 1995, A&A, 294, 792. 
Hollenbach, D.J. & Tielens, A.G.G.M. 1997, ARA&A, 35, 179. 
Hollenbach, D., Kaufman, M.J., Bergin, E.A., Melnick, G.J. 2009, ApJ, 690, 
1497. 

Mai-coni, A. Testi, L. Natta, A. & Walmsley, C. M. 1998, A&A, 330, 696. 
Melnick, G. J., Genzel, R., Lugten, J. B. 1987, ApJ, 321, 530. 
Noel, B., JobHn, C, Maillard, J. R & Paumard, T. 2005, A&A, 436, 569. 
Le Petit, R, Nehme, C., Le Bourlot, J. & Rouelf, E. 2006, ApJS, 164, 506. 
Lis, D. C. Serabyn, E., Keene, J., et al. 1998, ApJ, 509, 299. 
Lis, D. C. & Schilke, P 2003, ApJ, 597, L145. 
Offer, A.R. & van Dishoeck, E.R 1992, MNRAS, 257, 377. 
Offer, A.R., van Hemert, M.C. & van Dishoeck, E.E 1994, J.ChemPhys 
100, 362. 

Parmar, PS., Lacy, J.H. & Achtermann, J.M. 1991, ApJ, 372, L25. 
Pilbratt, G. L. et al. 2010, A&A, 518, LI. 
Poglitsch, A, Waelkens, C, Geis, N. et al. 2010, A&A, 518, L2. 
Schoier, E. L. van der Tak, E. F. S., van Dishoeck, E. E. & Black, J. H. 2005, 
A&A, 432, 369. 

Simon, R., Stutzki, J., Sternberg, A. & Winnewisser, G. 1997, A&A, 327, L9. 
Storzer, H., Stutzki, J. & Sternberg, A. 1995, A&A, 296, L9. 
Storey, J. W. V., Watson, D. M. & Townes, C. H. 1981, ApJ, 244, L27 
Tappe, A.. Lada, C. J., Black, J. H. & Muench, A. A. 2008, ApJ, 680, LI 17. 
Tielens, A.G.G.M. & Hollenbach, D. 1985, ApJ, 291, 722. 
van der Werf, P, Stutzki, J., Sternberg, A. & Krabbe, A. 1996, A&A, 313, 633. 
van der Wiel, M.H.D. et al. 2009, A&A, 498, 161. 
van Harrevelt, R., & van Hemert, M. C. 2000, J. Chem. Phys., 1 12, 5787 
Walmsley, C. M., Natta, A., Oliva, E. & Testi, L. 2000, A&A, 364, 301. 
Wampfler, S.E, Herczeg, G.J., Braderer, S. et al. 2010, A&A, 521, L36. 
Young Owl, R.C., Meixner, M.M., Wolfire, M., Tielens, A.G.G.M., Tauber, J. 
2000, ApJ, 540, 886. 



Goicoechea et al.: OH emission towards the Orion BarPDR 



Appendix A: analysis and intensity extraction 

The best-fit OH radiative transfer model was obtained by finding 
the minimum ;i'^-value defined as 

where n is the number of observed OH lines, p is the number 
of free parameters in the fit, /' , and /' , are the observed and 

^ ' ODS mod 

modelled line integrated intensities, and cr|^^^ is the Icr uncer- 
tainty of /^^^. Therefore, we fit absolute intensity values and not 
line intensity ratios. Figure IaTI shows the best grid of radiative 
transfer models with a density of rin-lO^ cnT^ and a beam fill- 
ing factor 77=0.1, defined asrj - q-^^, where Qj and ^Ipsp are 
the sky solid angles subtended by the source and by the PACS 
PSF respectively. Figure lATI shows the best source-averaged OH 
column densities, where N{Oli)heam~tl x N(OH)source- 

The list of OH line intensities around the position where the 
higher excitation OH lines peak is shown in Table IA.II In order 
to better sample the spectrometer PSF in these extended emis- 
sion observations, the OH line intensities were computed by 
co-adding the measured fluxes over four adjacent spaxels (the 
[3,4], [3,3], [2,4] and [2,3]). Owing to the broader PSF width 
at long wavelengths, the entire array was used to extract the 
OH ~163 yum line intensities. 




5 10'* 10'° 1.5 10'° 2 10'^ 



Fig. A.l. xijxlod a function of T^ and A^(OH) for a grid 
of models with rin-lO^ cm"^ and 77=0.1. The set of parameters 
giving the minimum values with a confidence level of 99.9% 
are shown. 



Table A.l. PACS OH line intensities towards the Orion Bai" 
PDR''. 





OH 


lobs (C obs) 


(p.m) 


transition 


(xlQ-^ erg s-'cm-^sr"') 


119.441 


^Hj/z y=5/2+-3/2- 


10.09 (0.21) 


119.234 


2n3/2 y=5/2--3/2+ 


9.44 (0.21) 


84.597 


^Hj/z y=7/2--5/2+ 


4.07 (0.79) 


84.420 


2n3/2 y =7/2+ -5/2- 


2.38 (0.91) 


163.396 


^Hi/z y=3/2--l/2+ 


1.34 (0.17) 


163.015 


^Hi/z y=3/2+-l/2- 


2.74 (0.16) 


79.179 


'ni/2-2n3/2 y=i/2+-3/2- 


9.17(1.30) 


79.115 


'ni/2-'n3/2 y=i/2--3/2+ 


9.65 (1.54) 


65.279 


2n3/2 y=9/2+-7/2- 


1.02 (0.49) 


65.131 


^ny2 y=9/2--7/2+ 


1.32 (0.54) 


Near the OH emission peak: a20oo-5*35 


"'21.9',52ooo^-5°25'06.7 
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Appendix B: Online figures 




Fig.B.l. Rotational energy diagram of '^OH showing the Hnes 
detected by HerschelfPACS (in microns). Detected transitions 
are shown in black while other pumping transitions discussed in 
the text are shown in red. OH ^113/2 and rotational ladders 
are produced by the spin-orbit interaction, while the A-doubling 
splitting of each rotational level is produced by the nuclei rota- 
tion and the unpaired electron motion. The A-doubling splitting 
has been enhanced for clarity. Hyperfine structure is not shown. 





Fig.B.2. Same as Figure |2]for the OH ^01/2-^03/2 7=1/2^3/2 
cross-ladder lines at ~79 fim, the OH ^03/2 J-5/2^3/2 ground- 
state lines at ~119/^m, and the OH 'ni/2 7=3/2— »l/2 excited 
lines at ~163/im. Note the decrease of OH line intensity with 
distance from the ionization front. 



7 



